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We present a QCD sum rule calculation of the spin-symmetry violating universal 
function Xi{ v • v'), which appears at order 1/mg in the heavy quark expansion of 
meson form factors. This function vanishes in the standard approximation, where 
radiative effects are neglected. For the first time, the complete set of diagrams arising 
at order a s is evaluated. In particular, we find = —(3.8 ± 0.4)% at zero recoil, 

indicating that 1/mg corrections induced by the chromo-magnetic moment operator 
are small. 
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I. INTRODUCTION 



Recently, the discovery of a spin- flavor symmetry that QCD reveals for heavy quarks 
has increased the prospects both for a reliable determination of some of the weak mixing 
angles and a study of nonperturbative QCD in semileptonic decays of heavy mesons 
These symmetries are responsible for restrictive relations among weak decay amplitudes and 
reduce the number of independent form factors. The heavy quark effective theory provides 
a convenient framework in which to analyze such processes It allows a systematic 

expansion of decay amplitudes in powers of 1/m.Q. The coefficients in this expansion are 
universal functions of the kinematic variable y — v ■ v', where v and v' denote the velocities 
of the initial and final mesons. They are independent of the masses and spins of the heavy 
quarks, but characterize the properties of the light degrees of freedom. At leading order, 
only a single function £(y) appears. The conservation of the vector current implies that this 
celebrated Isgur-Wise form factor is normalized at zero recoil, allowing model-independent 
predictions unaffected by hadronic uncertainties. 

Already at order l/rng, however, one encounters a set of four additional universal func- 
tions £3(2/) and Xi(y) f° r i = 1)2,3, as well as a parameter A, which denotes the mass 
difference between a heavy meson and the heavy quark that it contains ]7). These functions 
break the spin-flavor symmetries and induce corrections to the relations valid in the infinite 
quark mass limit. At zero recoil, two of them can be shown to vanish, Xi(l) = X3(l) = 
0, but £3(1) and ^2(1) are nonzero and lead to observable effects. For instance, the leading 
correction to the ratio of the axial form factors describing B — > D*iv decays is given by || 



Am B m D , A 2 (gi ax ) _ 1 _ _A_ 
(m B + m D *) 2 Ai(g^ ax ) 2m c 



;i + 3r)&(l)+4(l-r)x2(l) , (1) 



where r = m c /mb- An understanding of the universal form factors is at the heart of nonper- 
turbative QCD, but is necessary for any phenomenological application of the heavy quark 
expansion. QCD sum rules are particularly suited for this purpose and have recently been 
employed to calculate the Isgur-Wise form factor [PHT2] , the mass parameter A [13||b§ , and 



the universal functions that appear at order I/tuq in the heavy quark expansion ||. For 
these subleading form factors an interesting result was obtained. Within the approximations 
usually adopted in QCD sum rules, one finds the parameter-free predictions 

6,(1) = 5zi = i, X2(y) = 0, (2) 

where D is the space-time dimension. Corrections to these relations could only come from 
radiative effects or from condensates of rather high dimension, both of which are expected 
to be small and are usually neglected. Provided those corrections are indeed negligible, the 
predictions (Q) allow for an essentially model-independent description of the 1/toq correc- 
tions at zero recoil. For instance, using r ~ 1/3 and A ~ 0.5 GeV, one finds ~ 0.9 for the 
right-hand side of (0). 

In this paper we investigate the corrections to the second relation in fl2|), which arise if 
one goes beyond the standard approximations by including diagrams of order a s in the sum 
rule for x^iv)- In particular, we calculate for the first time the two- loop radiative corrections 
to the perturbative triangle diagram. 
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II. DERIVATION OF THE SUM RULE 



The heavy quark effective theory provides an expansion of hadronic matrix elements in 
powers of I/tuq. Its construction is based on the observation that, in the limit ttlq ^> Aqcd, 
the velocity v of a heavy quark Q is conserved with respect to soft processes In this 
limit it is possible to remove the mg-dependent piece of the heavy quark momentum Pq by 
the field redefinition 

h Q (v,x) = e m ^ x P + Q(x) (3) 

where P + = + is a positive-energy projection operator. For simplicity we will abbrevi- 
ate h = hg(v,x). The new field h carries the residual "off-shell" momentum k = Pq — m,QV, 
which is of order Aqcd- The effective Lagrangian for the strong interactions of the heavy 
quark becomes PJI1,|T5 



C mET = hiv-Dh 



h(iD) 2 h + Z^ha af3 G al3 h .... (4) 



where Z is a renormalization factor, D = d — ig s t a A a is the gauge-covariant derivative, and 
G is the gluon field strength. 

A similar expansion can be performed for a generic heavy quark current Q' T Q, with T 
being a 4 x 4 Dirac matrix carrying any number of Lorentz indices (T = 7 M (1 — 75) for the 
weak flavor- changing current). At tree- level, it simply reads 

Q'TQ ^Ch'Th + (5) 

where h' = fiQ/(v'), and C is again a renormalization factor ||16|]. 



A convenient way to evaluate hadronic matrix elements in the effective theory is by using 
the covariant trace formalism introduced in Refs. P,[l7|. Matrix elements of the leading 
current operator in (^) can be parameterized as (y = v ■ v') 

(M'{v')\h'Th\M{v)) = -£{y)tY{M'{v')rM{v)}, (6) 

where £(y) is the Isgur-Wise function, and 

v p ( — 75 ; pseudoscalar meson ,^ 

[V) y/ m + I ^ • vector meson ^ ' 

is a spin wave function associated with the pseudoscalar or vector meson M, which has 
the right transformation properties under Lorentz boosts and heavy quark spin rotations. 
At subleading order in the 1/rn.g expansion, matrix elements receive contributions from 
higher-dimension operators in the effective Lagrangian (|j) and in the effective current (|5|). 
A particular type of correction comes from an insertion of the chromo-magnetic operator 
contained in £hqet- The corresponding matrix elements of 

0{x) = i jdyT{[h'Th} x ,[Sfh a aP G aP h ] J (8) 

can be parameterized by a tensor form factor, 
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(M'(v')\0\M(v)) = -Atr{ X aP (v, v') M'(v')T P + ia a pM(v)}. (9) 

Because v a P + a a pM.{v) = 0, the most general decomposition of x a<3 consistent with Lorentz 
invariance and the heavy quark symmetries involves two real, scalar functions \2 and x% 
defined by 

x <#( Vl v >) = (r, V - «'V) X2(y) - 2ia^ xM- (10) 

Since we have factored out A in (Q), the functions Xi{v) are dimensionless. 

The QCD sum rule analysis for these subleading form factors proceeds along the same 



lines as that for the Isgur-Wise function. The procedure is outlined in detail in Refs. |8,|TT,H 
and we shall adopt the same notations here. We are interested in the analytic properties of 
the three-current correlator 

'dx dy e W-*- k -v) (0 | T{ [ qT M ,ti L,O(0),[hT M q] y }\0) 

= E 2 (u, u', y) tr{(i; V - u'V) T M ,P' + T P + ia a/3 P + T M } 

+ ^{u,u\y)ti{2a^T M ,P' + Y P + a aP P + T M }, (11) 

where P' + = |(1 Depending on the choice T M = —75 or T M = 7^ — v^, the heavy-light 

currents interpolate pseudoscalar or vector mesons, respectively. In the effective theory, the 
coefficients Hj are analytic functions in the "off-shell energies" uo = 2v-k and u' = 2v'-k', with 
discontinuities for positive values of these variables. They receive a double-pole contribution 
from the ground-state mesons M and M' associated with the heavy-light currents. Using 
the fact that the external momenta are P = itlqv + k and P' = itlq/v' + k', one finds 
that the pole positions P 2 = m 2 M and P' 2 = m 2 M , correspond to uo = to' = 2A, where 
A = rriM — rriQ = tum 1 — rnQ' 0- The residues are proportional to the universal functions 
Xi{y) an d Xz{y)i respectively. It follows that [|] 

Sr to ( W) a/,y) = ? 9 ^y A 9U ., i < = 2,3, (12) 

(oj — 2A + is) (lu' — 2A + ie) 

where F is defined by ( | qT h \M) = |Ftr{T Ai}. It is the analog of the meson decay 
constant in the effective theory |TT| . 

For large negative values of u and uj', the correlator can be calculated in perturbation 
theory using the Feynman rules of the heavy quark effective theory ||. The idea of QCD 
sum rules is that, at the transition from the perturbative to the nonperturbative region, 
nonperturbative effects can be accounted for by including the leading power corrections 
in the operator product expansion of the three-point function. They are proportional to 
vacuum expectation values of local quark-gluon operators, the so-called condensates |L8| . 
One then writes the theoretical expression for Hj in terms of a double dispersion integral, 

H* h (u;, uj' y) = fdudu'- — -r - — — - + subtractions, (13) 

J [y — uj — ie) \y' — uj 1 — ie) 

and performs a Borel transformation 
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n — > oo 
— — > oo 



with respect to c<j and u/. This yields an exponential damping factor in the dispersion 
integral and eliminates possible subtraction terms. Because of the flavor symmetry it is 
natural to set the Borel parameters equal, r = r' = 2T. It is then convenient to introduce 
new variables ct»± = \{v ± z/) to obtain 



/" dcu+ e" w+/T /" dcu_ 2pf 1 (cu++w_, y). (15) 



Following Refs. fll| , |12|| we perform the integral over U- and employ quark-hadron duality 
to equate the integral over uj + up to a threshold u Q to the Borel transform of the pole 
contribution in (|T2"D. This gives the QCD sum rules 



UJQ 

X i(y) F 2 A e- 2K ' T = J dco + e—l T p*{u> + , y) ; • = 2, 3. (16) 
o 

The spectral densities p* h arise after integration of the double spectral densities over 

In this paper we concentrate on the sum rule for Xi{y)- It has been shown in Ref. || 
that within the standard approximations, where one includes the bare quark loop as well 
as the leading nonperturbative corrections proportional to the quark condensate and the 
mixed quark-gluon condensate, there is no contribution to this function. The leading terms, 
then, are of order a s and come from the two-loop radiative corrections to the quark loop, 
the one-loop radiative corrections to the quark condensate, and the gluon condensate. In 
general, the calculation of all diagrams arising at order a s in the sum rule for a three-current 
correlator is rather tedious. In fact, the two- loop corrections to the triangle diagram have 
never been computed before. Because of the particular structure of the trace associated with 
S 2 in fln|), however, the calculation is considerably simplified in this case. Only the three 
diagrams shown in Fig. [I] contribute. Note that for each contribution the dependence of the 
spectral density p^ l {uj +) y) on u + is known on dimensional grounds, i.e. 

—pert 3 
P2 OC W + , 

oc const., 

pf G) cx % + ). (17) 

Instead of the spectral densities, it thus suffices to calculate directly the Borel transform 
of the individual contributions to E^ 11 , corresponding to the limit Uq — ■> oo in (|16"D. The 
uiq- dependence can then be deduced from flTT|). 

Let us now present the results of our calculation. The contributions from condensates are 
obtained by evaluating the diagrams shown in Fig. |T](b) and (c), using the one-loop tensor 
integrals given in Ref. M. We find 



econd (T ) = ^MI 

2 V ,yj 6tt 



r(y) 1 



y- i y + 1 



<S^(-U, (18) 
96tt \y + lJ v ; 
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where the function 



r(y) 



ln^ + vV-l) 



(19) 



satisfies r(l) = 1 and r'(l) = — |, such that there is no singularity in (|T£|) as y — > 1. The 
calculation of the two-loop diagram depicted in Fig. [l|(a) is more tedious. In an intermediate 
step, one encounters the two-loop tensor integral 



d D s d D t 



(s - t) a Sp t n 



(2tt) d (2tt) d (oj + 2v-s)(lu + 2vt)(uj> + 2v'-t) s 2 t 2 (s - t)< 



(20) 



For its evaluation it is convenient to represent the light quark propagators by Fourier trans- 
forms in a .D- dimensional Euclidean space (see |19[]) and use an exponential integral repre- 
sentation of the heavy quark propagators. After performing the integrations over the loop 
momenta, one is left with five parameter integrals, two of which become trivial after Borel 
transformation. The remaining three integrals stay finite for D = 4 and can be calculated 
in closed form. Our result is 



E^(T,y) 



a s T 4 
8tt 3 



y + 1 



1 - r(y) 

y-i 



(21) 



Based on ([Tq ) and (|T7|), we can now introduce back the dependence of the various contribu- 
tions on the continuum threshold ujq to obtain the final sum rule 



X2 (,)F 2 Ae- 2 ^ = -^(4- 

87r d \y + 1 



where 



a, 



T 



6n 



2 r 



1 - r(y) 
r(y) 1 



y-i 



y 



(a s GG) ( 2 



967r \y 



(t/ + i)' 



(22) 



5 n (x) 



1 

T(n + 1) 



dz z n e 



(23) 



III. EVALUATION OF THE SUM RULE 



For the QCD parameters in 



we take the standard values 

a s /-n = 0.1, 
(qq) = -(0.23 GeV) 3 , 
(a s GG) = 0.038 GeV 4 . 



(24) 



is 



The value of the strong coupling corresponds to the scale fi = 2A ~ 1 GeV, which 
appropriate for evaluating radiative corrections in the effective theory ]TT|. In addition, the 
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sum rule depends on the continuum threshold ujq and on the Borel parameter T. These 
quantities have recently been determined from the analysis of a QCD sum rule for the 
correlator of two heavy- light currents in the effective theory ||14|| . One finds good stability 



for ujq = 2.0 ± 0.3 GeV, and the consistency of the theoretical calculation requires that the 
Borel parameter be in the range 0.6 < T < 1.0 GeV. 

With this set of parameters one finds F = 0.30±0.05 GeV 3/2 and A = 0.50±0.07GeV for 
the hadronic parameters appearing on the left-hand side of (|22|). For reasons of consistency, 
however, one should rather work directly with the QCD sum rule for these quantities. This 



procedure minimizes the T-dependence of the final result. We thus use [14 



F 2 Ae -2A/T = ^! W^0\ _ ^o(gg) 



(25) 



where ~ 0.8 GeV 2 is the ratio of the mixed quark-gluon condensate and the quark 
condensate. Combining (|2~2"D and (p5|) , we express X2(y) as a function of u , T, and the 
QCD parameters. The numerical analysis of the resulting sum rule is shown in Fig. |2|. We 
observe excellent stability against variation of the Borel parameter T, and it supports the 
self-consistency of the approach that the stability region coincides with that determined in 
Ref. |T4il . Also the dependence on the threshold parameter is weak. Varying t^o between 1.7 
and 2.3 GeV changes X2{y) by no more than ±0.4%. For instance, we obtain from Fig. ^|b 
the zero recoil prediction %2(1) = —(3.8 ± 0.4)%. 

Finally, we note that a rough estimate of the form factor can be obtained by neglecting 
the nonperturbative terms in ( |22"D and (^). Independently of uq and T, this gives the simple 
function 



pert/ ^ _ a s ( 2 \ 2 [1 -r(y) 



xT\y) 



97T \y + 1 



(26) 



which is smaller than the complete sum rule result by roughly a factor |. At zero recoil, for 
instance, one finds xT* (1) — — ^("^s/ 71 ") — —2.6%. 



In summary, we have presented the complete order-a s QCD sum rule analysis of the 
subleading Isgur-Wise function X2{y)- We find that this form factor is small, typically of the 
order of —2% to —4%. When combined with the analysis of Ref. H, which predicts similarly 
small values for Xs{y)i but much larger values for the remaining two subleading form factors 
Xi(y) and £3(2/), these results strongly indicate that power corrections in the heavy quark 
expansion which are related to the chromo-magnetic moment operator are small. 
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FIGURES 



FIG. 1. Diagrams contributing at order a s to the sum rule for the universal form factor 
Xi{ v ' v ') : ( a ) perturbative contribution, (b) quark condensate contribution, (c) gluon condensate 
contribution. Heavy quark propagators are drawn as double lines, and a square represents the 
chromo-magnetic moment operator. 

FIG. 2. Numerical evaluation of the sum rules fl22] ) and (p5|): (a) form factor X2(v-v') evaluated 
for luq = 2GeV and 0.6 < T < 1.0 GeV; (b) dependence of X2(l) on the Borel parameter T for 
loo = 2.3, 2.0, 1.7 GeV (top to bottom). 
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